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Abstract 

We present the results of our calculations to a one, two, and three loop approx- 
imation of the e + e~— >e + e~ Bhabha scattering cross-section at small angles. All 
terms contributing to the radiatively corrected cross-section, within an accuracy of 
bo jo = 0.1%, are explicitely evaluated and presented in an analytic form. 0(a) and 
0(a 2 ) contributions are kept up to next-to- leading logarithmic accuracy, and 0(a 3 ) 
terms are taken into account to the leading logarithmic approximation. We define 
an experimentally measurable cross-section by integrating the calculated distribu- 
tions over a given range of final-state energies and angles. The cross-sections for 
exclusive channels as well as for the totally integrated distributions are also given. 
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"''Appeared as a contribution to the "Reports of the Working Group on Precision Calculations for the Z 
Resonance" CERN 95-03, March 1995. 
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1 Introduction 



An accurate verification of the Standard Model is one of the primary aims of LEP IJ. 
While electroweak radiative corrections to the s-channel annihilation process and to large 
angle Bhabha scattering allow a direct extraction of the Standard Model parameters, small 
angle Bhabha cross-section affects, as an overall normalization condition, all observable 
cross-sections and represents an equally unavoidable condition towards a precise determi- 
nation of the Standard Model parameters. The small angle Bhabha scattering process is 
used to measure the luminosity of electron positron colliders. At LEP an experimental 
accuracy on the luminosity of 

I— | < 0.001 (1) 
a 

will soon be reached 0. However, to obtain the total accuracy, a systematic theoretical 
error must also be added. This precision calls for an equally accurate theoretical expression 
for the Bhabha scattering cross-section in order to extract the Standard Model parameters 
from the observed distributions. An accurate determination of the small angle Bhabha 
cross-section and of the luminosity directly affects the determination of absolute cross- 
sections such as, for example, the determination of the invisible width and of the number 
of massless neutrino species N v ||. 

In recent years a considerable attention has been devoted to the Bhabha process 0, |J. 
However, the accuracy reached, is still inadequate. According to these evaluations the 
theoretical estimates are still incomplete, moreover, are in disagreement with each other 
up to 0.5%, far from the required theoretical and experimental accuracy 0. 

The process that will be considered in this work is that of Bhabha scattering when 
electrons and positrons are emitted at small angles with respect to the initial electron and 
positron directions. We have examined the radiative processes inclusively accompanying 
the main e + e~— >e + e _ reaction at high energies, when both the scattered electron and 
positron are tagged within the counter aperture. 

We assume that the center-of-mass energies are within the range of the LEP collider 
2e = y/s = 90 — 200 GeV and the scattering angles are within the range ~ 10 — 150 mrad. 
We assume that the charged particle detectors have the following polar angle cuts: 

0i < 0_ = Mr = 9 < 9 3 , 9 2 <9 + = ptf 2 , < 4 , 0.01 < 0, < 0.1 rad , (2) 

where p±, py, (p 2 , p%t ) represent the momenta of the initial and of the scattered electron 
(positron) in the center-of-mass frame. 

In this paper we present the results of our calculations of the Bhabha scattering cross- 
section with an accuracy of O(0.1%). The squared matrix elements of the various exclusive 
processes inclusively contributing to the e + e~— >e + e~ reaction are integrated in order to 
define an experimentally measurable cross-section according to suitable restrictions on the 
angles and energies of the detected particles. The various contributions to the electron and 
positron distributions, needed for the required accuracy, are presented using analytical 
expressions. 
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In order to define the angular range of interest and the implications on the required 
accuracy, let us first briefly discuss, in a general way, the angle-dependent corrections to 
the cross-section. 

We consider e + e~ scattering at angles as defined in Eq. (|2|). Within this region, if 
one expresses the cross-section by means of a series expansion in terms of angles, the 
main contribution to the cross- section da/d9 2 comes from the diagrams for the scattering 
amplitudes containing one exchanged photon in the t-channel. These diagrams, as it is 
well known, show a singularity of the type 9~ 4 for 9 — > 0, e.g. 

d9 2 

Let us now estimate the correction of order 9 2 to this contribution. If 

% ~ r«(i + c^) , (3) 

then, after integration over 9 2 in the angular range as Eq. (0), one obtains: 



~ C'(l + ci€ ; >-p). (4) 

min 

We see that, for 9 min = 50 mrad and max = 150 mrad (we have taken the case where the 
9 2 corrections are maximal), the relative contribution of the 9 2 terms is about 5xlO~ 3 Ci. 
Therefore, the terms of relative order 9 2 must only be kept in the Born cross-section 
where the coefficient c% is not small. In higher orders of the perturbative expansion 
the coefficient c\ contains at least one factor a/ir and therefore these terms can safely 
be omitted. This implies that, within our accuracy, only radiative corrections from the 
scattering type diagrams contribute. Furthermore, one should take into account only 
diagrams with one photon exchanged in the t-channel, since, according to the generalized 
eikonal representation, the large logarithmic terms from the diagrams with the multi- 
photon exchange are cancelled. 

Having, as a final goal for the experimental cross-section, the relative accuracy as in 
Eq. ([]]), and by taking into account that the minimal value of the squared momentum 
transfer Q 2 = 2e 2 (l — cos 9) in the region (|2]) is of the order of 1 GeV 2 , we may omit in the 
following also the terms appearing in the radiative corrections of the type m 2 /Q 2 with m 
equal to the electron (m e ), or the muon (m^) mass. 

The contents of this paper can be outlined as follows. In Section 2 we discuss the 
Born cross-section da B , taking the Z° boson exchange into account, and compute the 
corrections to it caused by the virtual and real soft photon emission. We present the 
results, as discussed above, in the form of an expansion in terms of the scattering angle 
9. We then define an experimentally measurable cross-section, <7 cxp , which is obtained by 
tagging the scattered electron and positron within a suitable range of polar angles and 
energies. We introduce the ratio £ = (x cxp /o"o by normalizing <j cxp with respect to the 
cross-section ctq — 47ra 2 / 'e 2 9\. In Section 2, by using a simplified version of the differential 
cross-section for the small angle scattering, we discuss the contribution to <r cxp from the 
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single bremsstrahlung process. In Section 3 we find all corrections of 0(a 2 ) to cr cxp caused 
by two virtual and real photon emissions. In Section 4 we consider 0(a 2 ) corrections 
caused by e + e~ pair emission. In Section 5 we derive the leading logarithmic corrections 
to the a 3 order by using the structure function method. In Section 6 we estimate the 
contributions of the neglected terms. Finally, in Section 7 we give the results obtained in 
terms of the ratio £ as functions of the experimental parameters. 

A more detailed derivation of these results will be presented elsewhere ||. 



2 Born and one-loop soft and virtual corrections 

The Born cross-section for Bhabha scattering within the Standard Model is well known 

da B a 2 



dQ 8s 



{AB 1 + (l-cYB 2 + (l + cYB 3 } , (5) 



where 
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Bi = (J) 2 |l + (^ 2 -^| 2 , B 2 = |1 + {gl - gl) X 

b, = ^|i + f + + ^) 2 (fe + x)i 2 + + f + - ^ a ) 2 (f e + x)T , 

As At 

X = s-m 2 + iM z T z ' ^ = t - Mf ' 

G M 2 1 1 

A = — ^ = (sin2^)" 2 , 9a = - 9v = - -(l-4sin 2 0j, 
2y lira 2 Z 

s = (pi+p 2 ) 2 = 4e 2 , t = -Q 2 = (p 1 - Pv ) 2 = ~s(l-c), 

c = cos 6 , 6 = pipy . 

Here 9 W is the Weinberg angle. In the small angle limit 

e 2 e A 9 6 

c = cos = 1 1 h ... • 

2 24 720 

Expanding the result (5) we have 

da B Ma 2 . 6 2 9 4 . 

where e = a/s/2 is the electron or positron initial energy. The contribution from weak 
interactions 5 woak , connected with diagrams with Z° boson exchange, is given by the 
expression: 

Q2 Qi 

5 woak = 2g 2 v f - — {gl + g 2 a ) Re X + — (g$ + g 4 a + Qg 2 v g 2 a ) \x\ 2 • (7) 

From Eq. (7) one can see that the contribution cf of the weak correction # weak into 
the coefficient C\ in Eq. (3) is 

w < 9 2 , {9l + 9l) M Z 2 (g^ + gl + Qglgl) M 2 

Ci ^ 2g v H — + fl max — j^- ~ 1 . (8) 
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According to our previous discussion after Eq. (4) this means that the contribution 
connected with the Z° boson exchange diagrams does not exceed 0.3%. We shall therefore 
neglect the radiative corrections to weak contributions, since they could contribute at most 
with terms < 10 -4 . 

In the pure QED case one-loop radiative corrections to the Bhabha cross-section were 
calculated a long time ago |J. Taking into account the contribution coming from the 
emission of soft photons with energy less than a given finite threshold Ae as well, one 
obtains for da^n / dc in one-loop approximation 

drr {1) da B 

^ = %^(l + ^ tual + U , (9) 
dc dc 

where da B ED is the Born cross-section in the pure QED case ( i.e. it is equal to da B with 
9a = 9v = ) and 

a 4e 2 6 e f sin2 1 dx 

cUtual + <Lf t = 2 - {2 [1 - In (— ) + 2 In ctg -] In — + / — ln(l - x) 

7r m l 2 Ae J cos 2 § x 

23 11 , Ae 2 ,, a 1 r vr 2 . 4 3 . 

1 In (—7)} + - ~, ^7 — (2c 4 - 3c 3 - 15c) 

9 6 m 2 ' f Ti (3 + c 2 ) 2 L 3 V ; 

Q 

+2 (2c 4 - 3c 3 + 9c 2 + 3c + 21) In 2 sin - 

6 6 
-4 (c 4 + c 2 - 2c) In 2 cos - - 4 (c 3 + 4c 2 + 5c + 6) In 2 tg- 

2 6 6 

+ - (11c 3 + 33c 2 + 21c +111) lnsin- + 2 (c 3 - 3c 2 + 7c - 5) In cos- 

O Zi Zi 

+2 (c 3 + 3c 2 + 3c + 9) St - 2 (c 3 + 3c)(l - c) S s ] . 

The value St (S s ) is defined by the contributions to the photon vacuum polarization func- 
tion n(t) [II(s)] as follows: 

H(t) = -(S t +h-l) + \ ft 2 L , (10) 

7T 3 9 4 7T 



ln^ , Q 2 = -t = 2e 2 (l-c) , (11) 



where 

Q_ 2 

m 2 

and we have taken into account the leading logarithmic part of the two-loop corrections 
to the polarization operator. In the Standard Model St contains contributions of muons, 
tau-leptons, W-bosons, and hadrons: 

S t = 51 + 8 T t + 8f + 8? , 5 s = 5 t (Q 2 -> -s) , (12) 

and the first three contributions are theoretically calculable: 

3 mi 9 



1 1 2\ 1 v r + 1 1 2 8 /, 4m 2 1 

r«i 1 / 2 A , ^W + l 1 9 11 /, 4m 2 

*r = r »«-4)ln— --^ + -, u. = (l + -^ 
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The contribution of the hadrons 5^ can be expressed through the experimentally mea- 
surable cross-section of the e + e~ annihilation ||. 

In the limit of small scattering angles we can present Eq. (9) in the following form: 

^ = ^£[i_n(tr 2 (i + 5), (14) 

5 = 2- [2(1 -£)lnl + -L-2] + - 6 2 A e + - 6 2 In A , 

7T A 2 7T 7T 

3 ,2 7 , 19 1 rN 

A e = —I 2 H / h - - 6.) , 

e 16 12 18 4 1 * s) ' 

. Ae . . Q 2 _ 2 
A = — , I = In — ~ In — . 
e s 4 

This representation gives us a possibility to verify explicitly that the terms of a rela- 
tive order 9 2 in the radiative corrections are small. Taking into account that the large 
contribution proportional to In A disappears when we add the cross-section for the hard 
emission, one can verify once more that such terms can be neglected. Therefore in higher 
orders we will omit the annihilation diagrams as well as multiple-photon exchange dia- 
grams in the scattering channel. The second simplification is justified by the generalized 
eikonal representation for the amplitudes at small scattering angles 
Let us introduce now the dimensionless quantity 



2 ri 

exp 

) 



4-7T or 

with Ql = e 2 6\ where <r cxp represents the experimentally observable cross-section: 



01 

47TCT 



J dx\ J dx 2 9(xix 2 - x c ) J d 2 qi 6{ J d 2 q 2 Q c 2 (15) 
dx\d 2 q\dx 2 d 2 q 2 

where and are, respectively, the energy fractions and the transverse components 
of the electron and positron momenta in the final state, sx c is the experimental cut-off 
on their squared invariant mass sx\x 2 , and the functions Q\ which take into account the 
angular cuts in Eq. (2) are defined as: 



0c = e{ o 3 _ M) 0(M _ 9l ) , e c 2 = e(e, - 0(M _ e 2 ) . (ie) 

x\t x±e ~ x 2 e x 2 e 

We restrict ourselves further to the symmetrical case only: 

e 2 = e 1 , e 4 = e 3 , p = e -±>i. (17) 

t>i 
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Let us define S as a sum of exclusive contributions: 



E = S + £ 7 + S 27 + S c+C + S 37 + S c+C 7 , (18) 

where S stands for a modified Born contribution, £ 7 , £ 27 , and S 37 stand for the con- 
tributions of one, two, and three photon emissions (both real and virtual), S c+C and 
S c+C 7 represent the emission of virtual or real (soft and hard) pairs without and with 
the accompaining real or virtual photon. 

By integrating Eq. (|6|) with the use of the full propagator for the t-channel photon, 
which takes into account the growth of the electric charge at small distances, we obtain: 



2 d6 2 



So = e\ J ^-(i - n(t))- 2 + e w + z e , (19) 

where 



9f 



2 

- I-) 2 I A 

J W — a l / ~~X7~ °wcak ) 



el 

is the correction due to the weak interactions and the term 



p 2 



comes from the expansion of the Born cross-section in Eq. (|5|) in powers of 9 2 . 

The one- loop contribution S 7 comes from one-photon emission (real and virtual). By 
adding to Eq. (14) the cross-section for hard emission calculated using a simplified version 



of the differential cross-section for small angle scattering [11] we obtain: 



p 2 , 

S 7 = £ / ^[1 - U(-zQl)]- 2 { J dx [{L z - l)P{x) (20) 

1 x c 



where 



P(x) = (\±^) + =\im{ 1 -±^9(l-x-A) + d + 2\nA)6(l-x)} , (21) 
1 — x a^o 1 — x 2 

is the non-singlet splitting kernel and 

k{x, z) = 2—4" I 1 + - *)\ + U + 6{x 2 p 2 - z) L 2 (22) 

+6{z-x 2 p 2 )L 3 . 
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Here L 7 = In ^—r- and 



Lt = In 



x 2 {z-\){p 2 -z) 
(x — z){xp 2 — z) 

L 3 = 



In 



(z — x 2 )(x 2 p 2 — z) 



111 



x 2 (x — z)(xp 2 — z) 
(z — x 2 )(xp 2 — z) i 
(x — z){x 2 p 2 — z) 



(23) 



3 Two-photon emission 

Let us now consider the corrections of the order of a 2 . They come from the two-photon 
emission as well as from the pair production, real and virtual. The virtual and soft real 
photon corrections can be obtained by using the results of Refs. [iT|-[n . 
The corresponding contributions to £ are: 



n S+V 



Eyy + Z VS + Z S S = (-) 2 / dz Z- 2 (l - U(-zQl))~ 2 R]7 

71 J 



(24) 



It is convenient to present the Rg + y m the following way: 



R 



5+V — r S+V + r S , +V(77) + r S+V-y 



= r s+vhl) = L 2 Z (2 In 2 A + 3 In A + + L z (-A In 2 A - 7 In A + 3£ 3 



rl +Vl = A[{L Z - 1) In A + -L z - l} 2 . 



(25) 



3 _45, 
2^ 2 16' 



Here the quantity A = 8e/e << 1 is the maximal energy fraction carried by a soft photon. 

The single hard photon radiation can be accompanied by real soft or virtual photons. 
It is useful to separate the cases of photons emitted by the same electron or positron or 
by both of them. 



d<r\ H ,s+v = da H{s+v) + d<r H (s+v) + dof s+v) + da ( x +v) 



(26) 



In the case where one of the fermions emits the hard real photon and another interacts 
with a virtual or soft real photon, we find: 



(S+V) _ " 



1-A 



I - 2 P ' 

2® [ dx—^- \ dzz- 2 \\-Tl{-zQ\ 

71 J 1 — X J 



2\l-2 



(27) 



{[l + e(x 2 p 2 -z)] (L z -l) + k(x,z)}(-)[(L z -l) In A + ^L 2 -1] . 

71 4 
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A more complicated expression arises when both photons interact with the same 
fermion. In this case the cross-section can be expressed in terms of the Compton tensor 
with a 'heavy photon' PT5| . We will consider below the case of the photon emission from 
the electron. An equal contribution arises from the positron activity. For the small angle 
detection of the final fermions we have R: 



^ HiS+V) = Z His+v) (28) 

l,a s9 fP 2 dz r 1 - A dx(l + x 2 ) r r .„ , A 3 X 

= -(-) / ^77 ; TvTvTn / — " -L z {(2 nA-nx+- 

2V ; ii ^ 2 (l-n(-zQf)) 2 ix c 1-x ziy 2 J 

[(L z - + 8) + k(x, z)] + - ln 2 x + (1 + 0) [-2 + In x - 2 In A] 

+ (1 - In x + 2 In A In x - In x ln(l - x) 

i9 r 1 -* dt , . . x(l - x) + 4x In x, (1-x) 2 , 

- n 2 x + / — In (1 - t) - 1 - '—} , 

Jo t y 1 2(1 + x 2 ) J 2(l+x 2 ) J ' 

where 9 = 8(x 2 p 2 — z). 

Let us consider now the contribution from the emission of two hard real photons. One 
can distinguish two cases: double photon bremsstrahlung, a) in opposite directions along 
electron and positron momenta, and b) in the same direction along electron or positron 
momenta. 

The differential cross-section in the first case can be obtained by using the factorization 



of cross-sections in the impact parameter space 15 . It takes the following form [11 



oo 1-A 1-A 

1 ,a s 



Sf = -(^.f J dzz- 2 [l-U(-zQl)]- 2 J dxi J dx 2 (29) 



x c £c 



1 — Xi 1 — X 2 

where 



' '' 1 ' ' 2 $(xi, z)${x 2 , z) 



$(x, z) = {L z - l)[8(z - l)0(p 2 -z)+ 9(z - x 2 )6(p 2 x 2 - z)\ 

+L 3 [-6(x 2 -z) + 9{z - x 2 p 2 )} + [L 2 + ^~ X l 2 ]e(z - x 2 )6(x 2 p 2 - z) (30) 

1 + x z 

+ - l)0i/ -,) + [«(!-,)- 0i, - n .n \ {^%-_% \ • 

Let us now turn to the double hard photon emission in the same direction, and the 
hard e + e~ pair production. We distinguish the cases of the collinear and semi-collinear 
kinematics of the final particles. In the first case all particles produced move in the cones 
within the polar angles 6$ « Q\ centred along the charged particle momenta (final or 
initial). The region corresponding to the case when both photons are radiated outside 
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these cones does not contain any large logarithmic contribution. In the semi-collinear 
region only one of the particles produced moves inside the cones, whilst the other moves 
outside them. 

In the totally inclusive cross-section the dependence on the auxiliary parameter 9$ 
disappears, and the total contribution has the form: 



Here 



E M = S HH = \{~f J dzz- 2 [l-n(-zQl)]- 2 J dx J dx x 



p 2 1-2A l-s-A 



I HH L Z 



X\{1 — X — Xi)(l — Xi) 2 

I HH = A 6{x 2 p 2 -z) + B + C 6[(1 - Xl ) 2 p 2 - z] 



(31) 



„,L Z , (p 2 x 2 — z) 2 . . 9 . , ,. , (1 — xi) 2 (l — x — xi) 

2 x 2 (p 2 x(l — Xi) — z) 2 xx\ 
B = ^ + ^ 2 {z - l)(p 2 -z)(z- x 2 ){z - (1 - x 1 ) 2 ) 2 (p 2 x(l - Xl ) - zf „ 



2 



2 >(p 2 x 2 -z)(z-(l-xi)) 2 (p 2 (l-xi) 2 -z) 2 (z-x(l-xi)) 

+ (x 2 + (l-xi) 4 )ln { }~ Xl?X \ +5 B , 
x(l — X — Xi) 

C = 1( 3(L Z + 2 ln |- 2 X( n (1 ~f l)2 v Z 2n ^ ^ " 2(1 " Xl)(3 ~ 2x(1 " Xlh 

(1 — Xi) z {p z X{l — Xi) — z){p z {l — Xi) — z) 

7 = l + (l-xi) 2 , (3 = x 2 + (1 - xi) 2 , 

J A — XXi(l — X — Xi) — X 2 (l — X — Xi) 2 — 2(1 — Xi)/9 , 

5b = xxi(l — x — xi) — x 2 (l — x — Xi) 2 — 2x(l — Xi)7. 
One can verify M that the combinations 



[ p2 ^£ r 77 , yH(S+V) , y,HH 

' 7r ! k z 2 [l-U(-zQl)f s+v + ^ 



a 2 /■/> ciz 
,a, n fp 2 dz 



V A z 2 [i-n(- 2 gfjp r ^ + E5+l/ + s ^ +s " ' (33) 

do not depend on A for A — > 0. 

The total expression S 27 which describes the contribution to Eq. (|i~8f ) from all (real 
and virtual) two-photon emissions is determined by the expressions in Eqs. (|24|), fl2~5|) , 
(0), (H), (H), and ©. 

Furthermore, it does not depend on the auxiliary parameter A and can be written as 
follows: 



S 27 = Xrr +V + 2 E^ +5 ) + Sf v+5) + S^ +5) + + 2£* H (34) 
= £ 77 + £ 7 + (-) 2 £<^ , £ = ln l^i . 
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The leading contributions E 77 , E 7 have the form: 



p i 

£77 = h?Lf f L 2 dz z -2 {1 _ n( _Q2 z )]-2 f dx {\p^( X ) [ 9{ X 2 ' f - z) + 1] 
2 7T J J 2 

1 

+ / jP(t)P£)0(ff?- z )} , (35) 



where 



p( 2 )(*) = f jP(t)P^)= lim{[(21nA + ^) 2 -46]5(l-a;) (36) 

1 + x 2 3 1 

+ 2 [ (21n(l - a;) -lna; + -) + -(1 + x)lnx - 1 + x] 6(1 - x - A)} , 

and 



oo 

1 ,OL 



2 z) 



E 7 = ity 2 J L 2 z dzz- 2 [l-Yl(-Q 2 lZ )Y 2 J dx x J dx 2 P(x 1 )P(x 2 )[9(z - l)9(p 

x c Xc/%1 

+6(z - x\)e(x\p 2 - z)\ [9(z - l)6(p 2 -z) + 6(z - x 2 2 )6(x 2 2 p 2 - z)] . (37) 

We see that the leading contributions to E 27 can be expressed in terms of kernels 
for the evolution equation for structure functions. The function 77 in the expression 
([34]) collects the next-to-leading contributions which cannot be obtained by the structure 
functions method |TJ], 13]. It has a form given explicitely in Ref. 0. 

4 Pair production 

In a similar way we consider also pair production. The corrections due to virtual e + e~ 
pairs can be extracted from Ref. |I2] . Using the expression for soft pair production cross- 



section [14| one obtains for the contribution of the virtual and soft real e + e pairs to E 



the following result: 



t+v = (^) 2 / ^' 2 [1 - U(-zQ 2 )]- 2 R^y , (3* 
i 



.2, . L , , 17 



Rs+v = In A + -) + L z (— - - In"' A 

2D 4 
-- In A --&) +0(1) 
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In this expression the quantity A = 5e/e <<lis the maximal energy fraction carried by 
a soft pair. 

Here we have taken into account only e + e~ pairs. An order of magnitude of the pair 
production correction is less than 0.5%. A rough estimate of the muon pair contribution 
gives less than 0.05% since hxQ 2 /m 2 ~ 31nQ 2 /m 2 . Contributions of pion and tau-lepton 
pairs give even smaller corrections. Therefore, within the 0.1% accuracy, we can omit any 
pair except e + e~ . 

Let us consider now the hard pair production. In this case we can restrict ourselves 
only to the collinear region where the produced pair moves in the small cones within 
the polar angles 9q around the fermion momenta. Indeed, the non-leading semi-collinear 
region will give contributions of one order of magnitude smaller than the leading ones and 
therefore within the required accuracy one can neglect them 0. 

The pair production contributions to £ appear from two regions with the pair com- 
ponents moving respectively along initial and scattered electrons and analogously for 
positrons. This separation allows us to carry out the integration over the energy frac- 
tion of the pair components. The resulting contribution to £ has, with both directions 
included, the form : 

P 2 1-A 

= -(^-f J dz z- 2 [l -U(-zQl)]- 2 J dx{R (x)[£ 2 [l + 6(x 2 p 2 -z)} (39) 

1 x c 

+ AC In x] + 2 9(xp - 1) 9{x 2 p 2 - z) C f{x) + 2/i(ar) C} , 

where we put the auxiliary parameter 9$ equal to 9 2 = z9f, and introduce the following 
notations: 



Ro(x) = \ ^^ + ^-^(4 + 7x + 4x 2 ) + 2(l + x)lnx, (40) 

O X X <oX 

131 136 2 2 4 20 2. „ r , 4 

fix) = 1 x x 2 ; + -(-4x z - 5x + 1 + - 

JK ' 9 9 3 3x 9(1 -x) 3 V x 

4 1 13 1 

H ) In (1- x) + -(8x 2 + 5x -7 ) lnx - 2 ln 2 x 

1 — x 3 1 — x 1 — x 

Ax 2 Y, In (1-3/) , 0/1 , _ , . W1 , , J , ln(l-y). 



/ dy [ V) +2(l + x)[6 + lngln(l-x)+ [ dy 1 y) Ul) 
J y J y 

o w o tf 



1 — x J y 
o y 

, , , 116 151 2 Ax 2 20 l r 2 5 

= -^r + ^r x + ^r + ^r-^n — t + tt8x 2 -iox-io + - — inx 

9 9 3a; 3 9(1 — x) 3 1 — x 

2 A A 

+ -[~Ax 2 - 5x + H 1 ] ln(l - x) - (1 + x) In 2 x 

3 x X x 

ln (! " V) , , , m 4 l l\ ln (! " 2/)/ 



/" , mil — y) , . 4 /■ , in (1 — y), , 

+2i + l ( 2 + /rfy^ ^ + lna;ln 1-x / rfy^ %2 

j V 1 — a; J y 

o o 

where f x {x) = -xf(~). 
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The total contribution of virtual, soft, and hard pairs does not depend on the auxiliary 
parameter A and can be written as follows: 



~ = + e£'" = ^) 2 £ 2 J dzz' 2 \\ - Tl(-zQl)]-* (43) 

1 

1 

/ dxR(x)[6(x 2 p 2 - z) + 1] + (-) 2 £0 C 

J IX 



k e e 

J 



where 



1 — x 2 

R(x) = 2 (1 + x) lnx + (4 + 7x + Ax 2 ) + -Pix) . (44) 

3x 3 

The explicit expression for C+C which collects the non-leading terms from Eqs. 
and (|39|) is : 



p 1 

C+C " =\j dzz~ 2 J dx[9(xp - l)9(x 2 p 2 -z)f + f 1 -^r- , (45) 

1 Xc 

where 

The same definition applies to f\. The regularization corresponding to the ( )_|_ 
prescription is as in Ref. [TO. 



5 0(a 3 ) corrections 

In order to evaluate the leading logarithmic contribution represented by terms of the type 
(aC) 3 we use the iteration of the master equation obtained in Refs. |L6|, [13| . 

To simplify the analytical expressions we adopt here a realistic assumption about the 
smallness of the threshold for the detection of the hard subprocess energy. By neglecting 
the terms of the order of: 



x c (-£) 3 ~ 10- 4 , (47) 

7T 

one should consider only the emission by the initial particles. Three photon (virtual and 
real) contribution to S have the form: 



13 



1 1 



£ 3 t = - (-£) 3 / dz z~ 2 [ dx! I dx 2 B(x x x 2 - x c ) [I 5(1 - x 2 ) P (3) (xi) (48) 

4 7T J J J 6 



1 X c X c 

6(x lP - l)6(x 2 lP 2 -z) + ^ 2 \x 1 )P(x 2 )9 l 6 2 ] [1 + 0(x 3 c )} , 
where P(x) and P^ 2 \x) are given by Eqs. ( plf ) and (|36| ) correspondingly, 

^2 = ^-4%^-.) , 

X 1 Jb 1 

p( 3 )( x ) = 5(1 - x ) A t + 0(1 - x - A) e t , 

A t = 48[^3-^2(lnA + ^) + i(lnA + ^) 3 ] , (49) 

# t = 48{- ] -^- f— - -f 2 + -ln(l -x) - -lnx + — ln 2 (l - x) 

* x 2 1 -x L 32 2^ 2 4 v ; 8 12 v ; 

+ — In 2 x — — lnxln(l — x)\ + -(1 + x) lnxln(l — x) — — (1 — 2) ln(l — x) 

X—x 

+ i-(5-3x)lnx-l(l-x)-i-(l+^)ln 2 x-i(l+x) / dy ^ 1 "^ } . 
32 16 32 8 7 y 







The contribution to X of the pair production accompanied by the photon emission 
when both, pair and photons, can be real and virtual is given below (with respect to Ref. 



IE] we include also the non-singlet mechanism of the pair production): 



p z 1 1 

1 ,OL 



£ c+c T=-(-£) 3 y dz z~ 2 J d Xl J dx 2 6{ Xl x 2 -x c ) 



1 Xc Xc 

{^[R P (x 1 ) - iiFfai)] 5(1 - x 2 )#(x 2 p 2 - z) + ~ P(x 2 )R( Xl ) 6M [1 + 0(x 3 )] , 

where 

i?(x) = R s (x) + -P(x) , iT(x) = ^—^(4 + 7x + 4x 2 ) + 2(1 + x) h\x , (50) 
3 3x 



i? p ( x ) = J R s (x)(- + 2ln(l-x)) + (l + x)(-ln 2 x-4 J dy- 



1 2 3 

+ -(-9 - 3x + 8x 2 ) In a: + -( 8 + 8x + 3x 2 ) . 

3 3 x 



y 



The total expression for E in Eq. p8| ) is the sum of the contributions given in Eqs. 
(|1^), fl20D , fl3"4|) , ([43]), fl48f) and d50|) . The quantity S is a function of the parameters x c ,p, 



2 



and Qi 
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6 Estimates of neglected terms 

Let us now estimate the terms not taken into account here in accordance with the required 
accuracy: 



a) Weak radiative corrections: 



aQl 
ttM? 



< 10" 5 . (51) 



b) Electromagnetic corrections to weak contributions, including interference terms : 



^w' ~ ^woakl^Oi— C < 10 4 . (52) 



7T 



Here 5 wcak is given by Eq. (7). 

c) Radiative corrections to the annihilation mechanism, including its interference with 
the scattering mechanism : 

rv 

EX ~ e\-C < KT 4 . (53) 

7T 

Our explicit expressions for E 7 , without annihilation terms, coincide numerically with 
the results obtained at the same order in Ref. |L7]| by using exact matrix elements. 



d) The interference between photon emissions by electron and positron: 



E int ~ 6\- < 10- 5 . (54) 

7T 



This contribution is connected with terms violating the eikonal form |T(| in the expression 



A(s } t) = A (s,t)e 1 ^ +0(— ) . (55) 

TCS 

e) The interference terms between one- and two-photon mechanisms of pair production, 
including the effect final particles identity: 

rv 

nr ~ (-) 2 ^ io" 5 . (56) 

7T 

f ) The semi-collinear mechanism of pair production gives a contribution which contains 
a small factor £ _1 with respect to the collinear terms, and is numerically small: Ej^ ir < 
10~ 4 . A more accurate estimate should be interesting. 
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g) The creation of heavy pairs (/i/z, rr, nn, ...) is at least one order of magnitude smaller 
than the corresponding contribution due to the light particle production and is therefore 
not essential. 

h) Higher-order corrections, including soft and collinear multi-photon contributions, 
can be safely neglected since they only give contributions of the type (a jC/ii) n for n > 4. 



7 Results and their discussion 

The total cross-section for the + c ~ distribution is: 

° = -qT e ' (57) 

where £ is given by Eq. ([181) . 

Let us define S° to be equal to S | n =o [see Eq. (19)], which corresponds to the Born 
cross-section obtained by switching out the vacuum polarization contribution IT(t) defined 
in Eq. (10). We obtain: 

a = £° (1 + S + F + 5 2 "< + 6 C+C ~ + 5 3 ^ + <5 C+C ~^) , (58) 



Qi 



where 



£° = £ |n=o = 1 " P + + S,| n =o - 1 - P' 2 , (59) 



and 



Zjq Zjq Zjq 

In the Tables below we give the values of 5p as well as their sum as functions of x c 
for the values of the parameters Q\ = 1.61°, and 82 = 2.80° defining the Narrow-Narrow 
(iVTV) case, and for 81 = 1.5°, and 8 2 = 3.15°, defining the Wide- Wide (W-W) case as in 
Ref. pj. 

Each of these contributions to £ has a sign which can be changed as a result of 
the interplay between real and virtual corrections. The cross-section corresponding to 
the Born diagrams for producing a real particle is always positive, whereas the sign of 
the radiative corrections depends on the order of perturbation theory for the virtual 
corrections: at odd orders it is negative, and at even orders it is positive. When the 
aperture of the counters is small the compensation between real and virtual corrections is 
not complete. In the limiting case of zero aperture only the virtual contributions remain 
giving a negative result. As a consequence we see that the radiative corrections for the 
iViV case are larger in absolute value than ones in the WW case. These corrections depend 
on 6*i, x c , and p. When x c — > 1 the corrections increase in absolute value. 
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Tables 



x c 


X 
Oo 


A7 



A 2 7 






A 3 7 






2^ 








NN 


p = 1.74 








0.2 


0.413E-1 


-0.922E-1 


0.506E-2 


-0.750E-3 


-0.561E-03 


0.911E-4 


-0.471E-1 


0.3 


0.413E-1 


-0.962E-1 


0.508E-2 


-0.633E-3 


-0.514E-03 


0.767E-4 


-0.509E-1 


0.4 


0.413E-1 


-0.101E+0 


0.452E-2 


-0.635E-3 


-0.483E-03 


0.669E-4 


-0.562E-1 


0.5 


0.413E-1 


-0.108E+0 


0.406E-2 


-0.672E-3 


-0.453E-03 


0.589E-4 


-0.637E-1 


0.6 


0.413E-1 


-0.119E+0 


0.385E-2 


-0.733E-3 


-0.414E-03 


0.517E-4 


-0.749E-1 


0.7 


0.413E-1 


-0.138E+0 


0.348E-2 


-0.823E-3 


-0.350E-03 


0.456E-4 


-0.943E-1 


0.8 


0.413E-1 


-0.174E+0 


0.397E-2 


-0.964E-3 


-0.249E-03 


0.442E-4 


-0.130E+0 


WW p = 2.10 


0.2 


0.409E-1 


-0.728E-1 


0.276E-2 


-0.579E-3 


-0.406E-03 


0.598E-4 


-0.301E-1 


0.3 


0.409E-1 


-0.771E-1 


0.290E-2 


-0.491E-3 


-0.377E-03 


0.506E-4 


-0.341E-1 


0.4 


0.409E-1 


-0.826E-1 


0.257E-2 


-0.502E-3 


-0.356E-03 


0.442E-4 


-0.399E-1 


0.5 


0.409E-1 


-0.902E-1 


0.223E-2 


-0.543E-3 


-0.333E-03 


0.387E-4 


-0.479E-1 


0.6 


0.409E-1 


-0.104E+0 


0.200E-2 


-0.606E-3 


-0.297E-03 


0.338E-4 


-0.620E-1 


0.7 


0.409E-1 


-0.127E+0 


0.207E-2 


-0.697E-3 


-0.238E-03 


0.304E-4 


-0.849E-1 


0.8 


0.409E-1 


-0.167E+0 


0.303E-2 


-0.838E-3 


-0.146E-03 


0.320E-4 


-0.124E+0 



Table 1: The values of 5 i in Eq. (||) for y/s = M z = 91.161 GeV, sin 2 9^ = 0.2283, 
Y % = 2.4857 GeV. The 'NN counter' corresponds to p = 1.74 and the 'WW counter' to 
p = 2.10 as a function of x c . 



x c 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


Rnn 


0.592 


0.627 


0.708 


0.770 


0.817 


0.763 


0.498 


Rww 


0.653 


0.740 


0.883 


1.020 


1.120 


0.994 


0.500 



Table 2: Values of Rnn and Rww where R represents the ratio of non- leading with 
respect to leading contributions and is defined as R = (- ) 2 ^4t- for y/s = M z = 91.161 
GeV, sin 2 6 w = 0.2283, T z = 2.4857 GeV. The 'NN counter' corresponds to p = 1.74 
and the 'WW counter' to p = 2.10 as above as a function of x c . 
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